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ABSTRACT 

Models that seek to unify BL Lacs and low-power radio galaxies predict that the 
two types of object should show similar isotropically emitted X-ray emission. Testing 
this is usually limited by difficulties in separating strong X-ray emission from a BL 
Lac nucleus and surrounding low-surface brightness emission. In this paper we report 
ROSAT HRI observations of the z = 0.055 BL Lac object PKS 0521-365. We are able 
to separate a luminous extended X-ray component from the bright nucleus. Using a 
new radio map, we show that it is unlikely that the extended emission is due to inverse- 
Compton scattering of photons from the active nucleus, and instead interpret it as 
thermal emission from dense, rapidly cooling gas. This is a more extreme environment 
than is found in typical FRI radio galaxies, and may pose a problem for unified models. 
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1 INTRODUCTION 

In unified models for radio sources, BL Lac objects are the 
low-power radio galaxies whose jets are the most highly 
beamed along our line of sight (e.g. Browne 1983; An- 
tonucci & Ulvestad 1985; Ulrich 1989; Urry & Padovani 
1995). Their rapid variability, apparent superluminal mo- 
tion, strong point-like emission in radio, optical and X-ray, 
and the detection of some sources in 7-rays, are all explained 
if we are seeing emission from a relativistic jet closely aligned 
with our line of sight. Low-power radio galaxies then repre- 
sent the 'parent population' of unaligned objects whose jets 
are less favourably aligned. These low-power radio galaxies 
are likely mostly to be Fanaroff & Riley (1974) class I (FRI) 
objects, but the parent population may also include some 
transitional objects and low-excitation FRII radio galaxies. 

An important test of this unified model is the degree 
to which the isotropic (alignment-independent) properties 
of BL Lacs are similar to those of the parent population 
of radio galaxies. Such tests have been made, on the whole 
successfully, by looking at the extended radio emission (e.g. 
Antonucci & Ulvestad 1985; Kollgaard et al. 1992; Perlman 
& Stocke 1993, 1994) and properties of the host galaxies (e.g. 
Ulrich 1989; Abraham, McHardy & Crawford 1991; Wurtz, 
Stocke & Yee 1996; Falomo 1996) although there is some 
evidence that there are too few BL Lacs associated with the 
most luminous host galaxies (Wurtz et al. 1996). 

Another isotropic indicator is the clustering environ- 
ment. Using two-point correlation analysis of optical fields, it 
has been shown that FRI radio galaxies are normally found 



in groups or clusters of galaxies (Longair & Seldner 1979; 
Prestage & Peacock 1988) and BL Lacs seem also to in- 
habit groups or poor clusters (Pesce, Falomo & Treves 1995; 
Smith, O'Dea & Baum 1995; Wurtz, Stocke & Ellingson 
1997) though it appears that, at least at low redshift, BL 
Lacs are not often found in the dominant galaxies of rich 
clusters (Prestage & Peacock 1988; Owen, Ledlow & Keel 
1996; Wurtz et al. 1997); for example, Owen et al. (1996) 
find no BL Lacs at the centres of Abell clusters, a result 
inconsistent at the 95% confidence level with the numbers 
expected from the unified models of Urry & Padovani. 

Clustering environment may also be investigated by X- 
ray observations. It has long been known that many objects 
close to the FRI-FRII luminosity boundary are associated 
with rich clusters having luminous X-ray haloes. Recent ob- 
servations with ROSAT have shown that more typical FRI 
radio galaxies have extended thermal X-ray emission on 
scales characteristic of groups or poor clusters (Worrall & 
Birkinshaw 1994). This offers a new way to test the unifi- 
cation hypothesis; such emission should be isotropic, and so 
we expect all BL Lacs to have X-ray haloes comparable to 
those of FRIs. This test is difficult because it requires us to 
separate any extended X-ray emission from the bright un- 
resolved emission of the BL Lac nucleus. In this paper we 
describe such an analysis of .ROSAT' observations of the BL 
Lac PKS 0521-365. 

PKS 0521—365 is a well-studied BL Lac with a redshift 
of 0.055, comparable to the redshifts of the radio galaxies 
studied by Worrall & Birkinshaw (1994). It is variously de- 
scribed in the literature as a blazar, a BL Lac object, or 
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an N-galaxy, and on multifrequency spectral index plots like 
those of Sambruna, Maraschi & Urry (1996) is placed among 
radio-selected BL Lacs. Its host galaxy is easily detectable 
in the optical [contributing ~ 50 per cent of the source lumi- 
nosity at 5500 A in an 8x8 arcsec effective aperture; Falomo, 
Scarpa & Bersanelli (1994)] and it exhibits strong, variable 
broad emission lines (Scarpa, Falomo & Pian 1995). The 
host galaxy is a giant elliptical (Wurtz et al. 1996). Pesce et 
al. (1995) suggest that the excess galaxy count around the 
object corresponds to a cluster of Abell richness or more; 
they identify at least one, and up to four associated galaxies 
within 100 kpc. However, the cross-correlation analysis of 
Wurtz et al. (1997) suggests a poorer cluster, with richness 
class < 0. 

In the radio, the source has a 408-MHz flux of 36.1 
Jy (Wright & Otrupcek 1990), corresponding to a power at 
that frequency of ~ 4 x 10 25 W Hz" 1 sr" 1 ; this places it 
slightly above the nominal FRI-FRII luminosity boundary 
(~ 1 x 10 25 W Hz -1 sr -1 at this frequency), though of course 
some of the 408-MHz emission is contributed by the core, 
presumed to be beamed. It exhibits a core-halo-hotspot mor- 
phology on arcsecond scales (Wardle, Moore & Angel 1984; 
Ekers et al. 1989; see also section ^), which, together with 
its comparatively high radio power, may suggest that it is 
an aligned version of a transitional FRI-FRII object. The 
prominent radio jet is also seen in optical synchrotron emis- 
sion, extending about 6 arcsec from the nucleus (e.g. Keel 
1986, Macchetto et al. 1991). No motion of the core com- 
ponents was detected in VLBI observations (Tingay et al. 
1996) and this, together with the comparatively low ratios of 
nuclear to host-galaxy optical emission (Falomo et al.) and 
radio core to extended radio flux (Antonucci & Ulvestad 
1985), suggests a source that is only moderately relativisti- 
cally boosted along the line of sight compared to the most 
extreme BL Lacs. It was for this reason that we selected 
it as a suitable candidate for an X-ray search for extended 
emission with the ROSAT HRI. PKS 0521-365 has already 
been extensively observed at X-ray wavebands, with Ein- 
stein (Worrall & Wilkes 1990), EXOSAT (Sambruna et al. 
1994) and the ROSAT PSPC (Pian et al. 1996), and was 
detected in 7-rays by EGRET (Lin et al. 1996), but none of 
these observations provides the resolution necessary to sep- 
arate compact from extended structure. The source appears 
to have a concave spectrum in the X-ray waveband; the en- 
ergy index is ~ 1.0 in the soft band, flattening to ~ 0.5 at 
higher energies. 

Throughout the paper we use a cosmology in which 
Ho = 50 km s" 1 Mpc" 1 , q = 0. At the distance of PKS 
0521—365 1 arcsec corresponds to 1.475 kpc. 



2 OBSERVATIONS 

We observed PKS 0521-365 with the ROSAT HRI for 46.6 
ks between 1995 Feb 09 and 1995 Feb 11. The data were 
analysed using the IRAF Post-Reduction Off-line Software 
(PROS). 

PKS 0521-365 is detected at 10620 ± 140 net counts 
in a 2.5-arcmin radius circle with background taken from 
an annulus between 2.5 and 3.5 arcmin. The source circle 
was chosen by expanding the radius until the background- 
subtracted counts enclosed stayed constant within the er- 



rors. A 25-arcsecond radius circle around a nearby point 
source (1 arcmin to the north) was excluded from the analy- 
sis. The detected counts correspond to a rate of 0.228±0.003 
counts s - . The source centroid is consistent to within a few 
arcseconds with the best available optical and radio posi- 
tions. 

A simple analysis, binning the data into 1-ks intervals, 
suggests that the source was not strongly variable during 
the 110 ks that contain our observations. To compare the 
source flux with the earlier ROSAT PSPC observations, we 
re-analysed the 4.8 ks of PSPC data from the ROSAT public 
archive, taking account of the nearby background sources. 
We detect the source at 3610 ±60 PSPC counts in the same 
source region, or 0.75 ± 0.01 counts s _1 . A spectral fit to a 
power- law model between 0.1 and 2.4 keV using the galactic 
N H (3.37 x 10 20 cm -2 ) of Elvis, Lockman & Wilkes (1989) 
gives an energy spectral index of 0.925 ±0.03 (x 2 = 23 with 
32 d.o.f.; errors are la for one interesting parameter), con- 
sistent with the results of Pian et al. (1996) and Sambruna 
(1997); the unabsorbed 1-keV flux density from the PSPC 
data is then 2.03 ± 0.03 /iJy, while the flux density from the 
HRI data on the same assumptions is 1.88 ±0.02 /ijy. If Nh 
is allowed to vary the best-fit spectral index is 1.10 ± 0.12 
with N H = i.QtoA x !0 2 ° cm" 2 ( X 2 = 17.8 with 31 d.o.f.; 
errors la for two interesting parameters). Incorporating the 
spectral uncertainties, the flux in the 0.1-2.4 keV ROSAT 
band seems to be 10 ± 2 per cent lower in the HRI data than 
it was at the epoch of the PSPC observations (1992 Aug 
29). We return to this point below. 



3 ANALYSIS 

Fig. [j] (left) shows that the source is approximately sym- 
metrical; radial profiling is therefore appropriate. To search 
for extension in PKS 0521—365 we first compared its ra- 
dial profile with the standard empirical expression (David 
et al. 1997) for the HRI point response function (PRF). It 
will be seen from Fig. ^ (left) that the radial profile is not 
well fitted by a standard ROSAT PRF (normalised to the 
central regions, which we expect to be dominated by the 
point source); there is a deficit of counts at small radii and 
an excess of counts at large ones. 

It is well known that images from the HRI can be 
badly affected by errors in the ROSAT aspect correction; 
the effect is to produce spurious elongation ('smearing') in 
the images of point sources, thus broadening the effective, 
azimuthally-averaged PRF. Such smearing is likely to be 
significant compared to the detector and instrument point 
response on scales of ^ 10 arcsec. As shown in Fig. [j], the 
source appears elongated, in its central regions, in position 
angle roughly 135°; aspect smearing could therefore be re- 
sponsible for some of the deviations from the nominal PRF 
shown in Fig. ^. 

The aspect errors are thought to arise because the star 
tracking system does not always calculate positions cor- 
rectly, due to variations in the gains of the pixels of the star 
tracker CCD. The spacecraft wobble exacerbates this prob- 
lem because it causes the guide stars to move across the 
tracking CCD on short timescales. Methods for correcting 
aspect smearing (e.g. Morse 1994) therefore rely on binning 
the data according to the phase of the wobble; so long as the 
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satellite roll angle and the properties of the CCD are con- 
stant, the aspect error should be a function of the wobble 
phase only. To do this, we applied iraf/PROS scripts pro- 
vided by J. Silverman and D. Harris, based on suggestions 
by G. Hasinger (for more details see Harris et al. 1998). We 
first established that all our data were taken at the same 
roll angle; this, along with the relatively small number of 
closely spaced observation intervals, allowed us to analyse 
all the data together. The scripts bin the data as a function 
of wobble phase, utilising the fact that the relative wobble 
phase at any time is given simply by the spacecraft clock 
(modulo clock resets, none of which occured during our ob- 
servation). We chose to divide the data into 20 phase bins, 
giving approximately 500 counts per binned observation; 
with 20 phase bins there are enough counts for centroid- 
ing in each bin, and examination of images derived from 
each bin showed them to be circular. A centroid was then 
found, using the detect suite of tasks within PROS, for each 
individual bin, and the images were restacked so that the 
centroids aligned. 

The resulting image (Fig. right) appears less ellipti- 
cal than the uncorrected image (Fig. nl, left) and is consid- 
erably better fit, in its inner regions, by the nominal HRI 
PRF (Figure || right); it appears that the stacking proce- 
dure has made the core slightly narrower than the nominal 
PRF, which is of course derived from observations that may 
themselves have been affected by aspect uncertainties. We 
conclude that the inner regions of PKS 0521—365 are well 
described as a point source. More importantly, even after 
dewobbling we still see an excess of counts over the PRF at 
radii of 10-40 arcsec. There is no reason to doubt that this 
corresponds to real extended emission. 

To characterise the scale size of this extended emission 
we fit models consisting of a /3-model (Sarazin 1986) and a 
point-source component, both convolved with the nominal 
PRF, to the background-subtracted radial profile of the de- 
wobbled image. To avoid large contributions to the fitting 
statistic due to the fact that the restacked data are narrower 
than the nominal PRF, we sampled the radial profile more 
coarsely in its central regions for these fits. We performed 
fits for various different values of the parameter (3 and the 
core radius. Figure |§] shows a plot of x 2 as a function of /3 
and core radius; the best fit (x 2 = 19 with 15 degrees of 
freedom) has /3 — 0.9 and a core radius of 8 arcsec, with 
a central normalisation of the /3-model of 5.4 ± 0.5 counts 
arcsec -2 . Figure ^ shows the best-fit model plotted with the 
data. Fitting a point source alone to the data gives unac- 
ceptably high values of the fitting statistic (x 2 = 124). 

For comparison, we fitted a similar range of /3-models to 
the archival PSPC dataset. Because of the broader intrin- 
sic PRF of the PSPC, smearing is not considered to be a 
problem with this instrument. The best-fitting point source 
and /3-model had a core radius of 35 arcsec with j3 = 0.9, 
as shown in Fig. ^. The data are better fit with this model 
than with a point-source model alone (best-fit x 2 = 3.6 with 
9 degrees of freedom), although the point-source fit is for- 
mally acceptable (x 2 = 6.9 with 10 degrees of freedom). A 
large range of possible core radii are acceptable, and the 
HRI best-fitting /3 model is among those consistent with the 
PSPC data within the 90 per cent confidence contour; the 
normalisations of the models are also consistent. If an appro- 
priately normalised /3-model with core radius 8 arcsec and 



— 0.9 is subtracted from the PSPC data the remaining 
emission is adequately modelled as a point source with no 
significant support for an additional extended component. 
We therefore conclude that the extended emission is well 
represented by the HRI best-fit model. The best-fit point- 
source count rates in this model in the PSPC and HRI data 
are consistent to within a few per cent, suggesting that the 
inconsistency between the total PSPC and HRI count rates 
(section ^) is a result of our not having taken into account 
the two spectral components corresponding to compact and 
extended emission. 



4 DISCUSSION 

The /3-model fits to the HRI data described above imply 
that the extended X-ray emission has a core radius of 12 kpc 
and contributes 860 ± 90 counts to the total in a 2.5-arcmin 
radius circle. The corresponding 0.2-1.9 keV luminosity, as- 
suming a Raymond-Smith model with kT — 1 keV and using 
galactic Nh, is approximately 8 x 10 35 W. This is consid- 
erably more luminous than the extended emission seen in 
most of the radio galaxies of Worrall & Birkinshaw (1994), 
although the FWHM of 14 kpc is comparable to some of 
the smaller objects. The central cooling time, on the same 
temperature assumption, is 6 x 10 s yr, which, if the emis- 
sion were thermal, would imply rapid cooling, with expected 
mass deposition rates of ~ 6 Mq yr~ . The radius at which 
the cooling time is comparable to the Hubble time would be 
~ 20 kpc, larger than the core radius, which means that a 
simple single-temperature model of the extended emission 
is not self-consistent. Cooling takes place for any reasonable 
choice of temperature of the X-ray gas. 

The extended X-ray component detected by the HRI is 
comparable in size to the extended radio emission in PKS 
0521—365 which might be thought to suggest an explanation 
in terms of inverse-Compton (IC) scattering by the energetic 
electrons in the radio lobes of photons from the cosmic mi- 
crowave background radiation (CMB) and from the BL Lac 
nucleus itself (cf. Brunetti, Setti & Comastri 1997). Fig. ^ 
shows the X-ray image superposed on a 1.4-GHz VLA radio 
map made from data supplied by G.G. Pooley. Since the flux 
on the shortest baselines of this radio observation (made at 
BnA configuration) is equal to the single-dish Parkes radio 
flux (Wright & Otrupcek 1990) of 16.3 Jy, we can be con- 
fident that the extended structure is well represented. The 
largest angular size (LAS) of the radio emission is about 50 
arcsec, whereas the FWHM of the extended X-ray compo- 
nent is 10 arcsec, and Fig. |^ shows that the radio emission 
is quite asymmetrically distributed around the nucleus, in 
contrast to the symmetry of the X-ray emission. However, 
we note that inverse-Compton emission would come prefer- 
entially from the parts of the radio lobes closest to the active 
nucleus, so this in itself does not make an inverse-Compton 
model implausible. Since we do not know the unprojected ge- 
ometry of the source, we model the extended radio-emitting 
component crudely as a uniform sphere centred around the 
nucleus. The flux density from the core is 2.7 Jy [consistent 
with the value of Antonucci & Ulvestad (1985)] leaving 13.6 
Jy in extended emission. Projection effects make it difficult 
to assess the true size (and therefore volume) of the lobes, 
but the flux expected from IC scattering from an object of a 
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given synchrotron (radio) flux is only weakly dependent on 
the object's size, and so we assume a radius for the sphere of 
half the LAS. Because the frequency of a scattered photon 
is raised by approximately the square of the Lorentz factor 
(7) of the scattering electron, and because the electron en- 
ergy spectrum is dominated by electrons with 7 ~ 100 [i.e. 
electrons near the assumed low-energy cutoff in the energy 
distribution; cf. Carilli et al. (1991)] the most significant con- 
tribution in the X-ray band from the IC process is due to 
far-IR photons. We accordingly model the spectrum of the 
BL Lac's nucleus in this region as a power law with spectral 
index 0.78 normalised to the measured 100-pim flux of Im- 
pey & Negebauer (1988). Pian et al. (1996) argue that the 
Doppler factor in the nucleus is of order unity, based on the 
relatively low radio core prominence discussed in section [j] 
and on models for the nuclear X-ray and 7-ray emission, and 
we adopt this value in order to infer a luminosity from the 
observed flux; we neglect the effects of possible circumnu- 
clear obscuration, so that our source model is effectively an 
isotropically radiating central source surrounded by a uni- 
form, spherically symmetrical electron distribution. Using a 
modified version of the code of Hardcastle, Birkinshaw & 
Worrall (1998a) we find that the IC process in this model 
cannot produce enough X-rays to account for the observed 
flux unless the lobes have magnetic field strength approx- 
imately an order of magnitude weaker than the value im- 
plied by equipartition of energy between magnetic field and 
relativistic electrons (for electron filling factor unity); this is 
relatively much weaker than fields that have been inferred in 
other sources from observations of X-ray emission attributed 
to inverse-Compton or synchrotron-self-Compton emission 
(Harris et al. 1994; Feigelson et al. 1995; Hardcastle et al. 
1998a; Brunetti et al. 1998). At equipartition, the X-rays 
from the IC process contribute at most 2 per cent of the 
observed flux, or less if the filling factor is lower or there is 
a significant contribution from relativistic protons. 

As a result of this simple analysis, we believe that it 
is unlikely that the extended X-rays seen around the source 
are due in large part to inverse-Compton emission. Instead, 
it seems more likely that they are indeed thermal emis- 
sion from a rapidly cooling central region, and that PKS 
0521—365 inhabits a cooling flow; this means that its envi- 
ronment is significantly different from those of the FRIs ob- 
served by Worrall & Birkinshaw (1994), which tended to lie 
in less dense environments with much longer cooling times. 
This model is qualitatively consistent with the bright ex- 
tended emission line region seen in PKS 0521—365 (Bois- 
son, Cayatte & Sol 1989) and with the strong polarization 
asymmetry seen in the radio observations. 

We attempt to fit the radial profile of the source with a 
cooling flow model based on that described by Hardcastle, 
Lawrence & Worrall (1998b). In this version of the model, 
the outer regions of the source are fit with an isothermal f3 
model while the temperature and density inside the cooling 
radius (r coo i) are power law functions of radius. The electron 
density follows 
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while the temperature is given by 
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r C ooi < r 

is the core radius of the /3 model, and dinner 
is the inner limiting radius of the cooling flow, allowing us 
to avoid infinities at r — 0. n e i, n e 2 and n e 3 are scale elec- 
tron densities and T c i, T e 2 and T e 3 scale temperatures; T e 3 
corresponds to the temperature of the non-cooling gas. It 
is clear that matching temperatures and densities allows us 
to write all the scale factors in terms of n e3 and T e3 . The 
results are only weakly dependent on ri nllcl - so long as it is 
small and we fix it at a value corresponding to 0.01 arcsec 
in what follows. 

The parameters a and 6 set the slope of the density and 
temperature distributions; for a realistic cooling model we 
expect a > and b < 0, so that density increases and tem- 
perature decreases with decreasing radius. The ideal gas law 
implies that pressure goes as r~ < - a+b ^ . If the atmosphere is re- 
quired to be close to hydrostatic equilibrium, then matching 
mass as a function of radius inside and outside the cooling 
zone leads to a constraint on r C ooi as a function of a and b: 
, 1/2 

'cool 



3/3 - (a + b) 

and if we assume some law for the radial dependence of 
pressure (P tx r _1 is consistent with observations of other 
cooling regions) then the model has only five free parame- 
ters: j3, r coro , a, T C 3 and n o3 , the last being a normalisation 
parameter that can be determined in the fit. In addition, 
such a model can only be physically realistic if the cooling 
time at r coo i is comparable to the system lifetime, or the 
Hubble time: 



30 x 10 b VkBT 10 
Tcooi = -. r k> 2 x 10 years 

ftp (/"cool ) 

(Sarazin 1986) where r coo i is in years, ksT in keV and n p in 
cm -3 . 

We fit a range of representative cooling models to the 
data. As before, j3 was chosen from a small number of pos- 
sible values (0.5, 0.667, 0.9) while r corc ranged from 1-1000 
arcsec. We tried values of 1, 1.5, 2, 5 and 10 keV for fcsT e 3. 
a was allowed to take the values 1.5, 1.75 or 2.0. No model 
consisting of a cooling flow alone was a good fit to the data. 
A number of models consisting of a cooling flow and a cen- 
tral point source were comparable in goodness of fit to the 
combination of /3-model and point source discussed above, 
but most of these were ruled out by the constraint on cooling 
time. Those that were not have small core radii, low exter- 
nal temperatures and steep temperature and density power 
laws: for example, the model with fcsT o3 = 1.5 keV, a = 2.0, 
P — 0.667, r corc = r coo i = 40 arcsec is an acceptable fit to 
the HRI radial profile (\ 2 = 25). In this model just under 10 
per cent of the total counts (1000 ± 100) are assigned to the 
cooling flow component. The nominal cooling time at the 
core radius is 2 x 10 10 years; the mass deposition rate is ap- 
proximately 40M© yr _1 and the implied densities at 20 kpc 
correspond roughly with those required for pressure equi- 
librium with the cold gas inferred from observations of the 
extended emission- line region (Boisson et al. 1989). Though 
the details of the model may not be correct, this shows that 
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the extended emission of PKS 0521—365 can plausibly be 
modelled as a cooling flow of this type. A strong nuclear 
source is still necessary; the flux of the point-like compo- 
nent in the cooling-flow model is reduced by only ~ 1 per 
cent compared to that derived from the simple /3-model fits 
of section |j| 



5 CONCLUSIONS 

We have carried out the first separation of nuclear and 
galaxy-scale extended X-ray emission in a BL Lac object 
and found evidence that PKS 0521—365 inhabits a dense 
and, presumably, rapidly cooling region of X-ray emission, a 
much more extreme environment than those found for 'typ- 
ical' FRI radio galaxies; we have shown that it can be mod- 
elled as a cooling flow in low-temperature cluster gas. If 
this result were extended to other BL Lac objects, it would 
cause difficulties for models that seek to unify FRIs and BL 
Lacs, and might imply some causal relationship between a 
dense and rapidly cooling atmosphere and the BL Lac phe- 
nomenon. However, it may be that PKS 0521—365, with its 
high power and intermediate radio structure, is not represen- 
tative of the BL Lac class. Further observations are planned 
both to verify the thermal nature of the galaxy-scale halo 
and to see whether PKS 0521—365 is unusual among BL 
Lacs in this respect. 
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Figure 1. Left: contour plot of the ROSAT HRI data for the inner regions of PKS 0521—365, smoothed with a Gaussian of a = 1 arcsec. 
Contours at 0.4 X (1, 2, 4, . . .) counts arcsec -2 ; the lowest contour corresponds approximately to the 3<r (99.87 per cent confidence) level 
for this convolving Gaussian (see Hardcastle, Worrall & Birkinshaw 1998c). Right: as left, but showing the data after the dewobbling 
procedure had been applied. 
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Figure 2. Radial profile of PKS 0521—365; counts arcsec as a function of distance from the source centre. The nominal on-axis PRF 
of the ROSAT HRI, normalised to the total counts within 5 arcsec of the centre, is plotted for comparison (solid line). Left: before 
dewobbling procedure. Right: after dewobbling procedure. 
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Figure 3. \ 2 as a function of core radius of the /? model for fits of the radial profile of PKS 0521—365 to a model consisting of the sum 
of a (3 model and point-like emission, both convolved with an HRI PRF. Points marked are the results of fits while the lines are natural 
cubic splines drawn through the points. 16 radial bins were used. The best fit (with \ 2 = 19, 14 d.o.f.) has (3 = 0.9, core radius 8 arcsec. 
Plusscs indicate fits with j3 = 0.35, stars j3 = 0.5, circles j3 = 0.67, Xs j3 = 0.9 and squares (i = 1.5. 
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Figure 4. The radial profile of PKS 052f— 365 plotted with the best-fit combination of /3 model and point source (dot-dash line). The 
contributions of the separate model components arc shown as dotted lines. The lower part of the plot shows the residual between data 
and model expressed in units of the fitting statistic x- A background level of 6.8 X 10 -2 count arcsec -2 has been subtracted. 
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Figure 5. \ 2 as a function of core radius of the ft model for fits of the PSPC radial profile of PKS 0521—365 to a model consisting of the 
sum of a (3 model and point-like emission. Points marked are the results of fits while the lines are natural cubic splines drawn through 
the points. 12 radial bins were used. The best fit (with x 2 = 3.6, 9 d.o.f.) has f3 = 0.9, core radius 35 arcsec. Plusscs indicate fits with 
(3 = 0.35, stars /3 = 0.5, circles (3 = 0.67 and Xs f3 = 0.9. 
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Figure 6. 1.4-GHz radio map of PKS 0521—365 overlaid on a smoothed X-ray image. The radio map was made from a 10-minutc 
snapshot observation with the VLA in BnA configuration, taken on 1989 Feb 14. The Gaussian restoring beam is 4.09 X 2.03 arcsec, as 
shown in the bottom left of the figure; contours arc at 8 X (1, \/2, 2, 2\/2, . . .) mjy beam -1 ). Note the hot spot to the SE and the jet 
extending NW from the core. The (dewobbled) X-ray image has been smoothed with a Gaussian with a = 4 arcsec; black corresponds 
to 3.8 counts arcsec -2 . 
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